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Abstract 



We investigate the relation between Cartan decompositions of the unitary group and 
discrete quantum symmetries. To every Cartan decomposition there corresponds a quantum 
symmetry which is the identity when applied twice. As an application, we describe a new 
and general method to obtain Cartan decompositions of the unitary group of evolutions of 
multipartite systems from Cartan decompositions on the single subsystems. The resulting 
decomposition, which we call of the odd-even type, contains, as a special case, the concurrence 
canonical decomposition (CCD) presented in [5], [6], [7] in the context of entanglement theory. 
The CCD is therefore extended from the case of a multipartite system of n qubits to the 
case where the component subsystems have arbitrary dimension. 
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1 Introduction 

Decompositions of Lie groups have been extensively used in control theory to design control 
algorithms for bilinear, right invariant, systems with state varying on a Lie group. Once it is 
known how to factorize a target final state Xf as the product 



Xi, ...,X r . In quantum information theory, a factorization of the type (1) can be interpreted as 
the implementation of a quantum logic operation with a sequence of elementary operations. In 
this case, the relevant Lie group is the Lie group of unitary matrices of dimensions n, U(n). In 
general, a decomposition of the unitary evolution operator of the form (1) is useful to determine 
several aspects of the dynamics of quantum systems including the degree of entanglement (see 
e.g. [16]), time optimality of the evolution [11] and constructive controllability (see e.g. [8], 
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then the task of controlling to Xf can be reduced to the (simpler) task of controlling to the factors 



[14])- 



complete and elegant results can be obtained, which also have important physical implications. 
Decompositions of the unitary group U (n) for large n exist and can be applied to the dynamical 
analysis of high dimensional quantum systems. However, the information obtained with this 
study is rarely as useful and of direct physical interpretation as in the low dimensional cases. 
For multipartite systems, this motivates the search for Lie group decompositions constructed in 
terms of decompositions on the single subsystems. We shall construct such type of decomposition 
in the present paper. 

The main motivation for the study presented here was given by the recent papers [5] [6] , [7] . 
In these papers, a decomposition of U(2 n ) called the Concurrence Canonical Decomposition was 
obtained for a quantum system of n two level systems (qubits). Such a decomposition has the 
above mentioned feature of being expressed in terms of elementary decompositions on the single 
qubit subsystems. It is related to time reversal symmetry and this raises the question of what 
in general the relation is between quantum mechanical symmetries and decompositions. As we 
shall see here, the answer to this fundamental question is instrumental in developing a general 
method to construct decompositions of multipartite systems from elementary decompositions of 
the single subsystems. We shall develop a decomposition which we call of the 'odd-even type' 
that contains the concurrence canonical decomposition as a special case. 

The paper is organized as follows. In Section 2 we review the basic definitions and results 
concerning discrete quantum symmetries and Cartan decompositions of the Lie algebra su{n) 
and therefore the Lie group SU(n). We shall stress the important result that, up to conjugacies, 
there are only three types of Cartan decompositions which are usually labeled as AI, All and 
AIII. In Section 3, we investigate the relation between Cartan decompositions and quantum 
symmetries and establish a one to one correspondence between Cartan decompositions and a 
subclass of symmetries which we call Cartan symmetries. To every Cartan decomposition of the 
Lie algebra u(n) and corresponding Cartan symmetry there corresponds a decomposition of the 
Jordan algebra of Hermitian matrices of dimension n, iu(n) equipped with the anticommutator 
operation. This is described in Section 4. This is also the crucial fact used to develop the general 
decomposition of the odd-even type for multipartite systems in arbitrary dimensions in Section 
5. This decomposition is a Cartan decomposition and, in Section 6, we show how to determine 
its type (AI or All). The Cartan decomposition also leads to a decomposition of the evolution 
of any quantum system into the product of an evolution with antisymmetric Hamiltonian and 
one with symmetric Hamiltonian with respect to a Cartan symmetry. This result is discussed 
in Remark 4.1. 
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2 Background material 

2.1 Discrete Symmetries in Quantum Mechanics 

Given a quantum system with underlying Hilbert space ft, a quantum mechanical symmetry is 
defined (see e.g. [9] Chapter 7, [13] Chapter 4) as a one to one and onto map 6 : H — > H such 
that physically indistinguishable states are also mapped into physically indistinguishable states 
i.e. for every \ip ><G H and <p\ G R , 

e(e < * 1 |^>) = e i * i e(|V'>), (2) 

for some ^£ E. Moreover 6 preserves the inner product of two states, namely if \a >:= @\a >, 
then, for any two states \a > and \(3 > 

| < a\0 > | = | < a\/3 > |. (3) 

In this definition, we omit for simplicity the consideration of selection rules and assume that all 
the states are physically realizable. A detailed discussion of this point can be found in [9]. 
According to Wigner's theorem [15], every such operation can be represented as 

9 = e**U, (4) 

where is a constant, physically irrelevant, real parameter, and U is either a unitary operator 
or an anti-unitary one. Recall that an anti-unitary operator U, \a >^ \a >:= U\a > is defined 
as satisfying 

< 0\a >=< f3\a >*, (5) 

U(a\a > +c 2 |/3 >) = c\U\a > +c* 2 U\(3 > . (6) 

Once a basis of the Hilbert space H is chosen, an anti-unitary operator U can always be 
written as 

U\a>= XK\a>, (7) 

where K is the operation which conjugates all the components of the vector \a > and X is 
unitary. 

A symmetry 0, whether unitary or anti- unitary, induces a transformation on the space of 
Hermitian operators A as 

A — > QAQ' 1 := 0(A). (8) 

It is in fact easily verified that 9(A) is a linear and Hermitian operator. Moreover the eigenvalues 
of 9(A) are the same as those of A and a set of orthonormal eigenvectors are given by 6 lay >, 
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where \atj > is an orthonormal basis of eigenvectors of A. It can be proved [9], [13], that, up to 
a phase factor, 9(A) := QAQ" 1 is the only choice that guarantees 

| < a\9(A)\(3 >= | < a\A\p > |. (9) 

Description of the symmetry is usually done by specifying how 6 acts on Hermitian oper- 
ators rather than how acts on states. This is because, Hermitian operators represent physical 
observables and therefore the action of 9 on observables is typically suggested by physical con- 
siderations. For example, the space translation symmetry has to be such that 

9(x)=x-a, (10) 

for some constant a where x is the position operator. As another example, the parity or space 
inversion symmetry is defined such that 

9(x) = -x. (11) 

On the other hand, specification of 9 on an irreducible set of observables uniquely determines 
up to a phase factor [9]. Recall that an irreducible set of observables {Aj} is defined such that 
if an observable B commutes with all of the {Aj}, then B is a multiple of the identity. 

An observable H is said to satisfy a symmetry or to be symmetric with respect to if 

9(H) = H, (12) 

or equivalently 

QH = HO. (13) 
It is said to be antisymmetric with respect to if 

9(H) = —H «-> QH = -HQ. (14) 

A special type of symmetry is the time reversal symmetry. In classical mechanics a time 
reversal symmetry changes a system into one which evolves with time reversal trajectories. This 
suggests to define a time reversal symmetry in quantum mechanics so that 6 acts on the position 
x and momentum operator p according to 

6(x) = x, (15) 

9(p) = -p. (16) 
This implies that the corresponding transforms momentum eigenvectors \p > as 

0|p>=|-p>. (17) 
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If the system under consideration has no spin degree of freedom, then x and p form an irreducible 
set of observables and therefore (15) and (16) uniquely specify the transformation 6 on the state. 
Moreover, from the definition of angular momentum L := x X p, we obtain 

9(L) = 9{x) x 6(p) = -L. (18) 

For a system with spin angular momentum S, we impose by definition, according to (18) 

0(S) = -S, (19) 

and x, p, S form an irreducible set of observables. If \m > is an eigenvector of the (spin) angular 
momentum corresponding to eigenvalue m, we have 

Q\m >= | - m > . (20) 

From these specifications, it is possible to obtain an explicit expression of the time reversal 
symmetry for a system of ./V particles with spin operators S\,...,Sn. It is given (in a basis of 
tensor products of the eigenstates of the z— component of the spin operators) by (see [9]) 

@ = e -jr( s i,y+ s, 2,y + ---+ s N,y)j(^ (21) 

where Sj y is the y component of the spin operator corresponding to the j—th particle, j = 
1, N, and K is the conjugation operator (same as in (7)). 

2.2 Cartan involutions and decompositions of su{n) 

We discuss next the Cartan decompositions for the Lie algebra su(n). What we say could be 
generalized to general semisimple Lie algebras. We refer to [10] for more details. 
A Cartan decomposition of su{n) is a vector space decomposition 

su{n) = K®V, (22) 

where the subspaces K, and V satisfy the commutation relations 

[/C,/C]C/C, (23) 

[K,P]CV, (24) 

[P,V]C]C. (25) 

In particular, notice that K, is a subalgebra of su(n). A Cartan decomposition of su(n) induces 
a factorization of the elements of the Lie groups SU (n) . Let us denote by e c the connected Lie 
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group associated to a generic Lie algebra C Then, given a Cartan decomposition (22), every 
element X in SU(n) can be written as 

X = KP, (26) 

where K G and P is the exponential of an element of V. Moreover if A is a maximal Abelian 
subalgebra of su(n), with AQV, then one can prove that 

\J Kee KKAK* = V. (27) 

This implies that one can write P in (26) as P = K\AK\ with ATi G and A G e" 4 . Therefore 
every element X in SU(n) can be written as 

X = K 1 AK 2 , (28) 

with K\,K2 G and ieeA This is often referred to as KAK decomposition. 

A Cartan involution of su(n) is a homomorphism : su{n) — > su(n) such that # 2 is equal to 
the identity on su(n). Associated to a Cartan decomposition (22) is a Cartan involution which 
is equal to the identity on JC and multiplies by —1 the elements of V, i.e. 

6{K) = K, VK G JC, (29) 

0(P) = -p, VP G V. (30) 

Therefore, given a Cartan decomposition, relations (29) and (30) determine a Cartan involution 
9. Viceversa given a Cartan involution 6, the +1 and —1 eigenspaces of 6 determine a Cartan 
decomposition. 

According to a theorem of Cartan [10], there exist only three types of Cartan decompositions 
for su(n) up to conjugacy. More specifically, given a Cartan decomposition (22) there exists a 
matrix H G SU{n) such that /C' := HICH*, V := HPH*, where JC' and V fall in one of the 
following cases labeled AI, All and AIII 3 . 

AI 

K! = so{n), V = so{n) L , (31) 

where so(n) is the Lie algebra of real skew-Hermitian matrices of dimension n, so(n)- 1 - is the 
vector space over the reals of purely imaginary skew-Hermitian matrices. The corresponding 
Cartan involution, which we denote by 6j, returns the complex conjugate of a matrix, i.e. 

6 I (A):=A. (32) 

3 In the following definitions and in the rest of the paper the inner product < A, B > in su(n) is defined as 
< A, B >:— Tr(AB*) and it is proportional to the Killing form (see e.g. [10]) 
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All 



K? = 8pQ, V = ap^)\ (33) 



where, we are assuming n even, and sp(^) is the Lie algebra of symplectic n x n matrices i.e. 
the Lie algebra of skew-Hermitian matrices A satisfying 



AJ + JA T = 0. 



(34) 



The matrix J is defined as 



The corresponding Cartan involution On is given by 




(35) 



9 H {A) := JAJ- 1 = -JAJ. 



(36) 



AIII. In this case Kl is the set of all the skew-Hermitian matrices A of the form 



A 



( 



R 



s 



) 



(37) 



where R G u(p), S G u{q), p, q > 0, p + q = n and Tr(R) + Tr(S') = 0. V' is equal to £ ' . The 
corresponding Cartan involution is given by 



Several authors have proposed Lie algebra decompositions for su(n) that, although special 
cases of the general Cartan decomposition, are of particular significance in some contexts. For 
example, Khaneja and Glaser [12] (see also [4] for the relation of this decomposition with Cartan 
decomposition) have factorized unitary evolutions in SU (2 n ), namely unitary evolution of n two 
level quantum systems (qubits), into local operations i.e. operations on only one qubit and two- 
qubits operations. This result has consequences both in the study of universality of quantum 
logic gates and in control theory. In the latter context, one would like to decompose the task 
of steering the evolution operator to a prescribed target into a sequence of steering problems to 
intermediate targets with a determined structure. 

Another decomposition which is of particular interest to us is the Concurrence Canonical 
Decomposition (CCD) of su(2 n ) which was studied in [5], [6], [7] in the context of entanglement 
and entanglement dynamics. In this decomposition, K! and V' are real span of tensor products, 
multiplied by i, of n 2 x 2 matrices chosen in the set {^2x2; o~x, o~y, o~z}, where cr Xj y :Z are the 
x, y, z Pauli matrices. In particular, K! is spanned by tensor products with an odd number of 



6 in (A) := I Ptq AIj 



(38) 



where the matrix I VA is defined as the block matrix I } 
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Pauli matrices and V' is spanned by tensor products with an even number of Pauli matrices. 
It was shown in [5], [6] that for n even this decomposition is a Cartan AI decomposition and 
for n odd is a Cartan All decomposition. One of the primary goals of the present paper is 
to extend the CCD to the case of multipartite systems of arbitrary dimensions. The CCD was 
also used in [1], [2] to characterize the input-output equivalent models of networks spin ^, in a 
problem motivated by parameter identification for spin Hamiltonians. Generalizations of these 
results for networks of spins of any value, in view of the results presented here, will be given in 
a forthcoming paper [3]. 

3 Relation between Cartan decompositions and symmetries 

The results of [5], [6] associate to the Concurrence Canonical Decomposition a time reversal 
symmetry. In particular, there is a relation between the involution 9 corresponding to the CCD 
and the time reversal symmetry O in (21) (with N = n the number of spin assumed all equal 
to ^). This relation is given by 

9{A) = OAQ-\ VAesu{2 n ), (39) 

where the right hand side needs to be interpreted as composition of operators. It is also easily 
seen, using only the fact that the time reversal symmetry is antiunitary and the general formula 
(7), that, if 6 is the time reversal symmetry on observables iA, we have 

9{iA) := SiAO' 1 = -iQAO' 1 = -i9(A). (40) 

This rises the question of whether there is in general a one to one correspondence between 
symmetries 0, 9 (8), and Cartan involutions 9 and therefore Cartan decompositions. Also the 
question arises on whether formula (see (40)) 

9(iA) = -iO(A),VA € u(n) (41) 

is always valid. We shall investigate these issues in this section. We shall see that only a 
particular class of symmetries, which we call Cartan symmetries give rise to Cartan involutions. 

Definition 3.1 A symmetry is called a Cartan symmetry if and only if O 2 is equal to the 
identity up to a phase factor. 

Cartan symmetries have the property that applied two times to any state return the physical 
state unchanged. For example the time reversal symmetry and the parity (11) are Cartan 
symmetries while the space translation symmetry (10) is not a Cartan symmetry. 
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Whether or not a symmetry is a Cartan symmetry can be verified once we have its repre- 
sentation in a given basis i.e. (cf. (7)) 

&\a>=XK\a>, (42) 

where X is unitary and K is the identity if is a unitary symmetry and is the conjugation of 
all the components of \a > if 6 is antiunitary. 6 is a Cartan symmetry if and only if 

XX = e^/nxn, (43) 

for some <j> € R in the antiunitary case and X 2 = e l ^I nxn for some 4> € R in the unitary 
case. This is clearly of the particular orthonormal basis chosen. If is antiunitary and T is a 
unitary transformation which transforms one orthonormal basis into another and XK describes 
the action of the symmetry in one basis then TXT*K describes the action of the symmetry in 
the new basis. It is easily seen that if X satisfies (43) so does TXT* and an analogous fact holds 
for unitary symmetries. 

Generalizing the approach in [5] [6] we now give the following definition. 

Definition 3.2 The transformation induced by a symmetry on su(n) is defined as 

9(A) := 0.40- 1 . (44) 

Notice this definition is analogous to the one of symmetries 9 on observables (8) which we repeat 
here with different notations: 

9(iA) := @iA@-\ VA e su{n). (45) 

In order to give an expression of the induced transformation in a given basis, we consider 
the antiunitary and the unitary case separately. In the antiunitary case, if K is the conjugation 
= XK, 0- 1 = X*K = KX*, 9 which gives 

9(A) = XAX*. (46) 

Analogously, one obtains 

9(A) = XAX*, (47) 

in the unitary case. 

Theorem 1 The transformation 9 on su(n) induced by a symmetry (@, 9) is a Cartan invo- 
lution if and only if (Q, 9) is a Cartan symmetry. Moreover, if (@, 9) is antiunitary, we have 
\/A G su(n), 

9(iA) := -i9(A). (48) 
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Moreover, if (@, 9) is unitary, we have \/A € su(n), 

9(iA) := iO(A). (49) 

Proof. It is easily verified that 9 defined in (46) or (47) is a homomorphism. Moreover, assume 
is a Cartan symmetry. Then we calculate (in the antiunitary case and analogously in the 
unitary case) 

6 2 {A) = X(XAX*)X* = XXAX*X* = A, (50) 

where in the last equality we have used the fact that is a Cartan symmetry. Therefore the 
associated 6 is a Cartan involution. 

Conversely consider a Cartan involution 6 on su(n), induced by a symmetry 0. Then we 
want to show that is a Cartan symmetry. 

Since 9 must be of the type AI, All or AIII, we must be able to write it as (32), (36) or 
(38) up to conjugacy. In particular there exists a unitary T such that (case AI) 

9(B) = TT*BTT* ,\/B € su(n), (51) 

or such that (case All) 

9(B) = TJ^BTJ-^^B e su(n), (52) 

or such that (case AIII) 

9(B)=TI p<q T*BTI p<q T* 1 (53) 
in the AIII case 4 We take in the cases AI, All and AIII given by (cf. (46) and (47)) 

= Tf*K, (54) 

= TJT*K, (55) 

and 

@ = TI p , q T*, (56) 

respectively. It is easily verified that these are all Cartan symmetries, i.e. XX = I with 
X = Tf*, X = TJT* and X 2 = I nxn with X = TI p>q T*. Moreover the choice is unique, up to a 
phase factor which does not change the property of the symmetry of being a Cartan symmetry, 
as the set of matrices su(n) is an irreducible set of skew-Hermitian operators. This concludes 
the proof of the theorem. □ 



4 In the case AI in appropriate coordinates the involution is equal to conjugation. If T is the matrix that makes 
the change of coordinates, every B G su(n) can be written as B — TAT* for a unique A in su(n) and therefore 
A = T'BT. Now 8i(B) = TAT* and replacing A = T*BT, one obtains (51). The other cases are analogous. 
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Remark 3.3 The theorem could have been stated in a somewhat stronger form. In fact, the 
proof shows not only that the symmetry corresponding to a Cartan involution is a Cartan 
symmetry by also that it exists and is unique up to a phase factor. Therefore there is a one to 
one correspondence given by (44) (45) between Cartan symmetries and Cartan involutions and 
therefore decompositions. 

Remark 3.4 It follows from the proof of the theorem that antiunitary Cartan symmetries 
correspond to Cartan involutions of the type AI and All while unitary ones give rise to Cartan 
involutions of type AIII. 

4 Dual structures of u(n) and iu(n); Commutation and anticom- 
mutation relations 

In this section, we study the dual structure of the Lie algebra u{n) of skew-Hermitian matrices 
and the Jordan algebra iu(n) of Hermitian matrices equipped with the anticommutator oper- 
ation. We shall see that to a Cartan decomposition of u{n) there correspond a decomposition 
of iu(n) which we also call 'Cartan' where the role of the subspaces are possibly reversed. This 
correspondence is crucial in the development of general decompositions for multipartite systems 
developed in the next section. The situation is somehow different if we consider decompositions 
of the type AI and All and if we consider decompositions of the type AIII. Therefore we shall 
consider the two cases separately. Only the case AI and All will in fact be used in the next 
section. 

Consider a Cartan decomposition of su(n) (22), (23), (24), (25) of the type AI or All, its 
corresponding Cartan involution 6 and Cartan symmetry 8 related through (48). This decom- 
position naturally extends to a decomposition of u(n) by replacing V with V © span{il nxn }. 
We shall denote this subspace, with some abuse of notation, again by V. So that 

u(n) = )CeV, (57) 

V = K,- 1 , where the orthogonal complement is now taken in u(n), and the commutation relation 
(23), (24) and (25) also holds, within u(n). The Cartan involution 6 of the type AI and All, is 
naturally extended to u(n) and span{il nxn } will belong to the —1 eigenspace of 6 so that the new 
definition of V is consistent with the fact that V is the —1 eigenspace of 0. The corresponding 
symmetry on iu(n) will be given by (48) or equivalently by (45). 

Consider now iu(n) 5 which has a structure of a Jordan algebra when equipped with the 
5 From now on in this section all the matrices are assumed skew-Hermitian so that matrices multiplied by i are 
Hermitian. An exception is in Remark 4.1 below where the H's denote Hermitian operators. 
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anti-commutator operation 

{iA, iB} := {iA)(iB) + (iB){iA). (58) 

Associated to a Cartan decomposition of u(n) (57) is a decomposition of iu(n), which we also 
call Cartan decomposition, given by 

iu(n) = iK® VP. (59) 

Moreover it follows from (45) that 9 is a homomorphism on the Jordan algebra iu(n). It is in 
fact an involution as 9 2 is equal to the identity map. It follows from (48) that VP and iK are, 
respectively, the +1 and —1 eigenspaces of 9 and therefore we have 

{VP, VP} C iV, {iP,iK} C iK, {ilC,i)C} C VP. (60) 

So the roles of the subspaces K and V are somehow reversed when going from u(n) to iu(n). 

In the case of a Cartan decomposition of the type AIII of su(n), the construction is similar. 
In this case, we extend the Cartan decomposition to u(n) by incorporating span{il nxn } into 
K rather than into P. The commutation relations (23) (24) and (25) are still valid with this 
modified definition. The induced decomposition on iu(n) given in (59) is such that iK and 
VP are the +1 and —1 eigenspaces of the involution 9. This follows from the correspondence 
between 9 and 9 which in this case is given by (49). We have 

{i)C,ilC} C iK, {iP,iK} C VP, {iP,iP} C iK. (61) 

Remark 4.1 (Decomposition of dynamics) It was pointed out in ([6]) that every evolution of 
a finite dimensional quantum system U := e lH can be decomposed as 

e iH = e m ae iH % ( 62 ) 

where the Hamiltonian H s is symmetric with respect to time reversal symmetry and the Hamil- 
tonian H a is antisymmetric (cf. (12)-(14)) with respect to time symmetry, i.e. 9(H S ) = H s and 
9{H a ) = —H a . Therefore every evolution can be decomposed into a time symmetric one and a 
time antisymmetric one. In view of the above treatment such a decomposition can be extended 
to any Cartan symmetry. For Cartan symmetries of the type AI and All, iH a and iH s are in 
the K and V subspace of the associated decomposition so that the decomposition (62) is Cartan 
decomposition (26). To this purpose also notice that the K Lie algebras in all the three types 
of Cartan decompositions correspond to semisimple compact Lie groups so that the exponential 
map is surjective [10]. The same argument can be repeated for Cartan symmetries of the type 
AIII with the only change that this time iH s <G K and iH a G P. 
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5 Cartan decompositions for multipartite systems in arbitrary 
dimensions; Decompositions of the odd-even type 



In this section we shall generalize the Concurrence Canonical Decomposition to the general case 
i.e. to the case of a multipartite system consisting of any number of quantum systems of any 
dimension. We shall call the general decomposition a decomposition of the even-odd type because 
the two subspaces in the Cartan decomposition consist of elements which are tensor products of 
an odd or even number of elements in appropriate subspaces. In doing this, we shall make use 
of the correspondence between decompositions in u(n) and decompositions in iu(n) described 
in the previous section. In particular, we shall consider decompositions of u(n)-iu(n) associated 
to antiunitary Cartan symmetries. To this correspond Cartan decompositions and involutions 
which have the property to extend to Cartan decompositions and involutions for multipartite 
systems as we shall now describe. 

Consider a multipartite quantum system composed of N quantum systems of dimensions 
n2,...,njv and with Hilbert spaces Hi,...,Hn- The space of skew-Hermitian (Hermitian) operators 
acting on the space TLj, j = 1, ...,N is u(rij) (iu(rij)). The space of skew-Hermitian (Hermitian) 
operators acting on the total Hilbert space Htot '■= 7~ti <8> 7~(-2 ® ■ ■ • <8> Hj\f is u{n\n2 ■ ■ ■ njv) 
(ra(niri2 • • • njy))- Consider now Cartan decompositions of u(nj), not necessarily all of the same 
type but all of the type AI or All, 

u(nj) = Kj ®Vj, (63) 

and the corresponding decompositions for iu(rij) 

iu(nj) = iJCj © iVj, (64) 

satisfying, with obvious modification of the notations, the commutation relations (23), (24), (25) 
and anticommutation relations (60). Let us denote by aj a generic element of an orthonormal 
basis in ilCj which is an Hermitian matrix. Also let us denote by Sj a generic element of an 
orthonormal basis in iVj which is also an Hermitian matrix. An orthonormal basis in iu(n\ri2 • 
■ ■ njv) is given by tensor products of the form 

F := Ti ® T 2 <g> • • • ® T N , (65) 

where Tj = o~j or Tj = Sj, with all the possible combinations of cr's and S"s in the places. 
We define I Q (l e ) the subspace of iu(niri2 • • - tin) spanned by tensor products which display an 
odd (even) number of elements a, so that we write 

iu{n\ni • ■ ■ njy) = Z Q © T e . (66) 
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We shall call this decomposition, along with the corresponding decomposition of u(n\ni • ■ • tin) 

u{n\ni ■ • • tin) = il © il e , (67) 
a decomposition of the odd-even type. We have the following result 

Theorem 2 The decomposition of the odd-even type (66) (67) is a Cartan decomposition which 
is associated to a antiunitary Cartan symmetry, i.e. 

[il ,il ] C il , [il ,il e ] C il e , [il e ,il e ] C il , (68) 

{1 ,2 } C l e , {l ,l e } C T , {J e , J e } C X e (69) 

Proof. The proof is by induction on the number of systems N. For N = 1, 1 D = iK,\ and l e = iV 
so that the commutation and anticommutation relations (68) and (69) are the same as (23)- (25) 
and (60), respectively. Assuming now (68) and (69) true for every number of subspaces strictly 
less then N we can verify (68) for N by using the formula 

[A®B,C®D] = ^{[A,C]®{B,D} + {A,C}®[B,D]), (70) 

and considering all the subcases. For example, to show the first one of (68) one considers the 
four cases, by indicating with the superscript the number of factors in the tensor products: 

CI: {iG 2^~\ B G ill, C € D G %L\ } 

C2: { A G T^-\ B G ill, C G Xf" 1 , D G il l a } 

C3: { A G Tf- 1 , B G ill, C G l?' 1 , D G iJ. 1 } 

C4: {i£ Tf- 1 , 5 G C G if" 1 , D G ^ } 

Analogously, one can verify (69) by using induction along with the formula 

{A ® B, C ® D} = ^([A, C] ® [B, D] + {A, C} © {B, D}). (71) 

□ 

Associated to a decomposition of the odd-even type is a Cartan involution on u(n\n2 ■ ■ ■ 
tin), 9 TOT , and the corresponding Cartan symmetry on the space iu(n\U2 ■ ■ ■ tin), 6 TOT . If 
0i, 6]\r and 9\, Ojy are the Cartan involutions and symmetries associated to the 1, 2, A^-th 
decomposition, Q TOT and Q TOT can be described as follows. 

Let A be an element of the orthonormal basis of il Q , i.e. it can be written as 

A = T 1 ®---®(iT k )®---®T N , (72) 
where Tj = Uj or Tj = Sj, with an odd number of cr's. Then 

6 TOT (A) = h{Ti) ® • • • ® k {iT k ) ® • • • <g> ^(Tjv) = ±A, (73) 
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since 9j(<jj) = -aj, 9j{Sj) = Sj, 9 k (ia k ) = ia k , and 9 k (iS k ) = -iS k . 

In general an element of the orthonormal basis in u{ri\n2 ■ ■ ■ un) is a tensor product of a 
and S elements, with i multiplying one of the elements. TOT is obtained by applying 9j in 
all the positions j without i and 9 k in the k— th position where there is the factor i. If in the 
k— th position there is a factor of the type a this gives a +lia factor when transformed. In the 
remaining terms, all the factors S are transformed into S while factors a give other factors of 
the type a and a collective factor (—l) p ~ 1 . Here p is the total number of cr's and this is 1 if p is 
odd and — 1 if p is even, so that TOT {A) = A in one case and 6 TOT {A) = —A in the other case. 

Analogously, one can treat the case where in the A;— th position there is a factor of the 
type iS. With a similar argument, one shows that Q TOT can be defined on tensor products by 
applying 9j in every j— th position which clearly gives a factor (— l) p where p is the number 
of factors a. This shows Q TOT and Q TOT are the involution and symmetry associated with the 
odd-even Cartan decomposition. 

An alternative treatment could have been to first define the involutions and symmetries and 
then to obtain the decomposition (66)-(69) in terms of eiegenspaces of these homomorphisms. 

Remark 5.1 The Concurrence Canonical Decomposition is obtained as a special case of the 
odd-even decomposition when the N systems are all two level systems and the Cartan decompo- 
sition chosen on each of them is of the type AIL This gives K, = su(2) = sp(l) and V = {0} in 
the decomposition of su{2) (22). It corresponds to a time reversal symmetry ((21) with N = 1 
and spin ^) which is indeed a Cartan symmetry. Notice that n = 2 is the only case where in the 
Cartan decompositions AI, All and AIII, we can take K, equal to the whole Lie algebra su(n). 
This fact makes it difficult, in higher dimensions, to obtain natural decompositions of dynamics 
into local and entangling parts as it was done for the 2-qubits case for example in [16]. 

6 The nature of the odd-even decomposition 

It follows (for instance) from formulas (68) (69) that the odd-even decomposition is a decom- 
position of the type AI or All, namely a decomposition corresponding to a Cartan symmetry. 
It is interesting to know how the choice of the single decompositions on the various subsystems 
determines whether the odd-even decomposition is of the type AI or AIL One reason for that 
is that one may want to further decompose the Lie algebra iX Q and therefore would like to know 
its nature. For example, it was shown in [6] that the Concurrence Canonical Decomposition is 
AI in the case of even number of qubit subsystems and All in the case of odd qubits. In our 
notation il is (conjugate to) 50(2^) for N even and sp(2 N ^ 1 ) for N odd. 
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In general this information can be obtained by a simple count of the dimensions. Recall that 
in a AI decomposition of u(n) the dimension of the Lie algebra /C in (57) is the dimension of 
so(n) i.e. di := n( - n ~^ while in a All decomposition of u(n) the dimension of the Lie algebra 
/C is the dimension of sp(^) i.e. dn := n ( n 2 +1 ) . These numbers are never the same and therefore 
they uniquely identify the type of decomposition obtained. We have the following result. 

Theorem 3 Consider an odd-even decomposition on N subsystem obtained by performing All 
decompositions on r subsystems and AI decompositions on N — r subsystems. Then the resulting 
decomposition is of type All if r is odd and of type AI if r is even. 

Proof. The proof is by induction on iV. If N is equal to 1 the result is obvious. Consider now 
N subsystems and consider first the case r odd. Assume, without loss of generality, that an 
All decomposition is performed on the iV— th subsystem. Let n\ (n 2 ) denote the dimension 
of the vector composed of the first N — 1 systems (of the N— th subsystem). By the inductive 
assumption the odd-even decomposition on the first N — 1 system is of the AI type. Therefore 
there are ni ( n 2 1-1 ) elements of the odd type i.e. tensor products containing an odd number of a 
matrices spanning the associated subalgebra /C. A remaining orthonormal set of n\ — " 1 (" 2 1 ~ 1 ) = 
ni(TH+i) even t yp g e i emen t s gp an t ne orthogonal complement in u(n\). The basis for the Lie 
algebra K, for the system composed of all the ./V subsystems is obtained by tensor products 
of the " l(n 2 1 " 1) odd elements with the n\ - " 2(w 2 2+1) = " 2(n 2 2-1) even type elements on the 
N— th subsystem or by products of the rai ( TO 2 1+1 ) e ven elements with the " 2< " 2 2+1 ) e ven type 
elements on the N— th subsystem. The dimension of the /C Lie algebra in the resulting odd-even 
decomposition is therefore 

ni(rai - 1) n 2 (n 2 - 1) | ni(ni + 1) ra 2 (n 2 + 1) _ nin 2 (nin 2 + 1) 

2 X 2 + 2 X 2 2 ' U4j 

which indicates a decomposition of the type All for the total system which has dimension nin 2 . 
An analogous reasoning proves the case where r is even. □ 

The result of [5] is obtained as a special case of the above theorem as in the case treated all the 
decompositions applied are of type All and therefore the Concurrence Canonical Decomposition 
is of type AI on an even number of subsystems and of type All on an odd number. 
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